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Abstract 

We evaluate, exactly in d, the master integrals contributing to massless 
three-loop QCD form factors. The calculation is based on a combi- 
nation of a method recently suggested by one of the authors (R.L.) 
with other techniques: sector decomposition implemented in FIESTA, 
the method of Mellin-Barnes representation, and the PSLQ algorithm. 
Using our results for the master integrals we obtain analytical expres- 
sions for two missing constants in the e-expansion of the two most com- 
plicated master integrals and present the form factors in a completely 
analytic form. 
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1 Introduction 



Recently the evaluation of the QCD form factors at the three-loop level has attracted 
much attention. The form factors constitute important building blocks for a number 
of physical applications. Among them are the two-jet cross section in e + e~ collisions, 
the Higgs-boson production in the gluon fusion and the lepton pair production in 
proton collisions via the Drell-Yan mechanism. The three-loop corrections to the 
form factors of the photon-quark and the effective gluon-Higgs boson vertex appear 
after integrating out the heavy top-quark loops. Let and Y^ v be the corresponding 
vertex functions. Then the form factors are defined by 

FM 2 ) = - 4( i- e)g 2 ^(^ r ^7M) , (1) 
F 9 (q ) = T g , (2) 

where q\ and — q 2 are the momenta of the incoming and outgoing particles (quarks, 
for the case of F g , and gluons, for the case of F g ), and q = q\ + q 2 is the momentum 
transfer. Here and below, if not stated otherwise, we put d = 4 — 2e. Within 
perturbative expansion, the form factors take the form 

* - i+E(£)'(fy>. m 



where Q 2 = —q 2 , and x is either g(quark) or ^ (gluon). One deals with the three-loop 
order and splits Fg into the singlet, fermionic and remaining gluonic part 

Ff = F^ + F^' nf + J2Q q 'Fi 3) ' sin9 , (4) 

where rif stands for the number of active quarks, and Q q is the charge of the quark 
q. The pole parts of Fq 9 and F^ in e were presented in Eqs. (3.7) of Ref. pQ 
and Eqs. (9) of Ref. |2j, respectively. The finite parts Fq 9 , F^ ,sm9 and Fg were 
presented in Ref. [3]. 

The integration-by-part reduction reduces the problem to the calculation of a 
small number of master integrals. All the master integrals apart from three most 
complicated master integrals contributing to the three-loop massless form factors have 
been evaluated in jH |5]. In fact, the word evaluated means here the evaluation up to 
the order of e which appears in the finite part of the form factors. Mathematically, 
this means the evaluation up to transcendentality weight six. About one year ago, 
one of the three most complicated master integrals (called A 9jl in jU [5J El [6] ) and the 
pole parts of Ag^ and Ag^ (shown in Figs. 1 and 2 in the next section) were evaluated 
analytically, while the e° parts of Ag t4 and Ag t2 were evaluated numerically — see 
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[31 [6|. Therefore, only the two (apparently, most complicated) pieces of the whole 
family of three-loop massless form factor master integrals are missing at the moment. 
Mathematically and aesthetically, it is desirable to obtain completely analytic results, 
and this is the problem we are going to solve in the present paper. 

Recently, in Ref. j7j a method of multiloop calculations based on the use of dimen- 
sional recurrence relations (DRR) [8] and analytic properties of Feynman integrals as 
functions of the parameter of dimensional regular izat ion, d, has been suggested. In 
the present paper we apply this method to evaluate, exactly in d, the master integrals 
contributing to massless three-loop QCD form factors. Using the derived expressions, 
we obtain analytic results for the missing two constants and thereby arrive at analytic 
expressions for the form factors. 

The key point of the approach of Ref. [7j is the analysis of the analytic properties 
of a given integral in a basic stripe of the complex plane d. The proper choice of the 
master integral, the basic stripe, and the summation factor can essentially simplify 
the analysis reducing the number of (or totally fixing) the constants parametrizing the 
homogeneous solution of DRR. The freedom of this choice, being an advantage, is also 
the only heuristic part of the method. For the case of massive tadpoles, this choice 
is relatively simple due to the possibility to get rid of the infrared and ultraviolet 
singularities by performing an analysis in an infrared-safe region d G (do, do + 2) and 
raising, if necessary, the powers of the massive denominators (see, e.g., example 2 
in Ref. [7]). For the case of massless on-shell vertex integrals this recipe does not 
necessarily work because raising the powers of the massless denominators also makes 
worse the infrared and collinear behavior of the integral. Thus, in this case, one 
should rely on an analysis of the corresponding parametric representation. A manual 
analysis of the parametric representation for the purpose of revealing the position and 
the order of the poles can still be a very complicated problem for the cases considered 
in this paper. Fortunately, the current version of the code FIESTA based on sector 
decompositions provides the possibility to solve this problem automatically. So, in 
order to apply the method of Ref. [7j to the calculation of a given master integral, 
we apply a complete set of various techniques: 

(i) a reduction to master integrals by two alternative ways: by a code based on 
[9] and the code called FIRE [10] to obtain DRR, 

(ii) a sector decomposition [11], [121 [13] implemented in the code FIESTA [131 EH] 
to determine the position and the order of the poles in the basic stripe, 

(Hi) the method of Mellin-Barnes representation [TSJ dH [T7] to fix the remaining 
constants parametrizing the homogeneous solution (if any), 

(iv) PSLQ [18] to guess the analytical expression for both the constants parametriz- 
ing the homogeneous solution and for the e-expansion of the master integral around 
d = A. 

As a result, we obtain representations for all master integrals in arbitrary d. The 
representations have the form of convergent series which allow, in particular, a fast 
high-precision calculation of the e-expansion around d — 4. 
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Figure 1: Master integrals for Ag^. 

The paper is organized as follows. In the next section we present an example of 
the calculation for the integral A-jp and give exact results for this integral and for the 
lower master integral We also present analytical expressions for the e-expansion 
of the integrals A 9 2 and Aj j4 . 

In the conclusion, starting from results for the form factors of Ref. [3] and sub- 
stituting the two constants by our analytic values, we present completely analytic 
expressions for the finite parts of the form factors. 
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Figure 2: Master integrals for Ag^- 

2 Master Integrals for Massless 
Three-Loop Form Factors 

Master integrals naturally form a partially ordered set. One master integral is said 
to be lower than the other master integral if the Feynman graph for the former 
can be obtained by contracting some internal lines from the Feynman graph of the 
latter. This ordering enables us to introduce the notion of complexity level of a 
given master integral which is the maximal number of nested lower master integrals. 
Owing to this definition, the master integrals with zero complexity level have no 
lower master integrals. The DRR for such integral is obviously homogeneous and its 
explicit solution is expressed in terms of T-functions. Moreover, it turns out that for 
three-loop on-shell massless vertex master integrals any integral expressed in terms 
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of T-functions has zero complexity level. We expect this situation to be general. 

Our primary goal is the calculation of the most complicated integrals, Ag j2 and 
Ag t 4 which are the last integrals in Figs. 1 and 2. However, in order to be able to apply 
the method of Ref. [7] we have to know all lower master integrals which are shown 
in the same figures. Four rows of diagrams in each figure correspond to complexity 
levels 0, 1, 2 and 3. Therefore, we start our calculation from the complexity level 1, 
then pass to the complexity level 2 and, finally, calculate the two master integrals of 
complexity level 3. Let us demonstrate an intermediate step of this procedure using 
the example of the integral A-j^. We directly follow the path of Ref. [7]: 

1. There are four lower master integrals, A4, A 51 , A 5 ^, and A 6 3 . Three of them 
are expressed in terms of T-functions, while the last one, A 6 3 , can be obtained 
using the same method with the final result conveniently represented as 

00 



A 6t3 (d) = A^(d)J2K' 2 3(d + 2k) + Al 3 (d), 

7T 4 2 11 - M CSc(^)cSc(f ) 



k=0 

Al'Ud) = -sm(nd)Al 3 (d) 



(3d-io)r(d-|)r(^i) 

>3 



.1,2^ (7d-18)sin(f)r(f-iy 

6 < sl ' ~ 37r»(d-3)r(f-3) ' [ ) 

2. Here and in what follows, we omit, for brevity, a power-like dependence of the 
master integrals on q 2 + iO which can easily be restored by power counting. 

Using the FIESTA program we determine the position and the order of the poles 
in the basic stripe which we choose as 5* = {d\ Red G (4, 6]}. The syntax for 
this analysis is SDAnalyze [U , F , h , degrees , order , dmin , dmax] , where U and F 
are the basic functions in the parametric integral corresponding to the given 
Feynman integral, h is the number of loops, degrees are the indices, order is 
the required order in e and dmin and dmax are values of the real part of d that 
determine the basic stripe. The output lists the values of d where the given 
Feynman integral can have poles. This feature appeared in the second version 
of FIESTA, but was not documented in [2] because testing was still in progress. 
So, after applying this procedure to A 7j2 we see3 that the integral has simple 
poles at d = 14/3, 5, 16/3, 6. 



4 In fact, the overall factor T(a~hd/2) where a is the sum of the indices is not taken into account by 
FIESTA but this can easily be done because the corresponding poles are explicit. Let us emphasize 
that FIESTA can report also on some fictitious poles. This can happen when contributions of 
individual sectors do have some additional poles which cancel in the sum. However, FIESTA itself 
can be used further to check whether the poles are indeed present or not. 
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3. The dimensional recurrence reads^l 

A 7 , 2 (d + 2) = c 7 , 2 {d)A 7t2 {d) (6) 
+c 6>3 (d)A 6>3 (d) + c 5>2 (d)A 5t2 (d) + c 5jl (d)A 5>1 (d) + c 4 (d)A 4 (d) 

where c n are some rational functions of d presented in the Appendix. 

4. Using the explicit form of the coefficient c 7j2 (d), we choose the summing factor 
as 

y _ (j - 3) COS (f) C°B (j - f ) COS (f + |) T (f - 9) 

= <<»= pjr^i ' (7 > 

Passing to the function A7 >2 (d) = 2(^)^2(0?), we obtain the following equation 

i 7 , 2 (d + 2) = A 7 , 2 (d) + A 6j3 (d) + A 5 , 2 (d) + A 5A (d) + A^d) , (8) 

where A n (d) = £(d + 2)c n (<f)A n ((i). The general solution can easily be con- 
structed using the explicit form of the integrals A 4 , A& t i, A 5)2 , and A 6i3 : 

00 

A 7j2 (d) = u{z) + J2 [^5,2(^-2-20 + i 5)1 (d _ 2 - 2/) + i* 3 (d-2-2Z) 
z=o 

00 00 00 

- J2 A H( d + 20 E i4g(d + 21 + 2fc) - ]T A 4 (d + 21), (9) 

Z=0 k=0 1=0 

where z = exp[i7rd]. 

5. The function £(d) has simple zeros at d = 14/3,5,16/3, therefore, A 7}2 (d) is 
regular everywhere in S except the point d — 6, where it has a simple pole. 
Besides, from the explicit form of the summing factor and from the parametric 
representation of A 7)2 (d) it is immediately clear that A 7)2 (d) grows slower than 
any positive (negative) power of \z\ when Imd — > —00 ( Imd — > +00). This 
fixes co(z) up to a function 

ai + a 2 cot(|(d-6)) ( 10 ) 

6. In order to fix the two remaining constants, we use data obtained from the 
Mellin-Barnes representation of A 7j2 (d) which can easily be obtained from the 
general Mellin-Barnes representation for the non-planar on-shell vertex diagram 



3 We use the integration measure d d k/ '{iir d / 2 ) per loop. 
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(see, e.g., Chap. 4 of [TT]): 




r (| - 2) r (| - 1) 2 v(d - 3)r(-z!)r(-z 2 )r(z 2 + 1) ; 



(27r) 2 7 J r(rf-2)r(f - 5) r(2d - 7)r(d - «i - 4) 



x 5 id Zl — rl r f^r - ^ - 6 ) r (^i + z * + i)r(d - zi - z 2 - 5) 



A 7)2 (d) 



r(f-, 1 -s) V2 

Xr - - Z 2 - r (-y + Zi + 2 2 + 7^ d*i d^ 2 . (11) 

Using the codes of Refs. [16J, at d = 6 — 2e and d = 5 — 2e we straightforwardly 
obtain 

A 7 , 2 (5-2e) = -_ + 0(e°). (12) 

Using these two values and also taking into account the fact that the singularities 
of the inhomogeneous part should be cancelled, we obtain 

7T 3 / TT 1Td\ TT 3 ( 7T 7T(A 7T 3 f ird IT 

uj(z) = = tan tan = tan 1 

V ; 20^ \10 2 J 36 \6 2 J 20^5 V 2 10 

TT 3 (ltd 7I"\ 7T 3 of Trd\ 137T 3 / Tld 

H tan 1 — H cot 3 — H cot — 

36 V 2 6/ 60 V 2 / 180 V 2 

7T 3 / 7r 7rcA n 3 f nd n\ , . 

H = cot = cot h - . (13) 

20^ V5 2 y 20^ V 2 5^ 1 ; 

Eqs. (jUJ), (TLB"]) , and (JTj) determine our final expression for v4 7)2 (<i). 

Two remarks are in order. First, our choice of the summing factor, the basic stripe 
and the master integral itself (we could have considered instead, e.g., an integral with 
some denominators squared and/or with numerators) may be not the most optimal 
one. With some other choice, we might have been able to fix the homogeneous part of 
the solution entirely within the method. However, given the number of the integrals to 
be considered and the absence of the general recipe for this choice, it was much more 
convenient to use in such cases additional data from Mellin-Barnes representations. 
In fact, for other integrals the number of the constants to be fixed was not greater 
than two. 

The second remark concerns the double sum in Eq. fl9]). Making a shift k — » k — l, 
we obtain the following triangle sum with the factorized summand: 

oo oo 

J2^:l(d + 2l)Y / AH(d + 2k). (14) 

1=0 k=l 
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The factorized form of the summand essentially simplifies the numerical calculation 
of the sum, making it possible to organize the calculations without nested do-loops. 
Proceeding in the same way for the rest of the integrals, we finally obtain general 
expressions for Ag^ and Ag t %. The resulting representations for arbitrary d are too 
lengthy to be presented here and can be obtained upon request from the authors. We 
present here only analytical results for the expansion of these two integrals around 
d = 4 which are most interesting for physical applications. The expansion for A) i4 
reads 



A 9A {4 - 2e) = e 



-37 B e 



+ 



+ 



+ 



1 8 

608C(3) 



+ 



1 + 



437T 2 

108" 



1 



+ 



109C(3) 14 53tt 2 
9 + ~9 + ~27~ 



9 



17 



311tt 2 481vr 4 



108 



12960 



949C(3) 2975vr 2 C(3) , 3463C(5) | g4 | Utt 2 | 



9 108 45 

434((3) 299tt 2 C(3) 3115C(3) 2 



9 

-339 



6 



+ 



77tt 2 



2539tt 4 
2592 



247613tt 6 
466560 



35vr 4 
108" 



7868C(5) 
15 



(15) 



For Asp, we arrive at the following result: 

2 5 



37b e 



A 9i2 (A - 2e) = e 

~31C(3) 



+ 



+ 



+ 



+ 



9e 6 6e 5 
50 181tt 2 
Y + 216 



+ 



20 17tt 2 



347C(3) 110 17tt 2 119tt 4 



9 



9 



432 



514C(3) 341tt 2 C(3) 2507C(5) 



9 36 
1516C(3) 737tt 2 C(3) 



9 



130 vr 

- + y 



24 

2 943vr 4 



1080 



+ 



15 

29C(3) 2 + 

195551tt ( 
544320 



170 19tt 2 163vr 4 



9 



6 



960 



2 , 2783C(5) 



6 

0(e) 



(16) 



3 Conclusion 



Eqs. (I15p and (fT6|) enable us to present completely analytic results for the three-loop 
corrections to the form factors defined by Eqs. (JTH4]) . Starting from Eqs. (8-10) of 
Ref. [3] and taking into account our analytic values of the e° terms in f)15p and (lT6"j) 
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we obtain the following analytic expressions 

61vr 2 C(3) 



->(3),g+n/ 



fin 

130017T 2 



72 

46594C(5) 



+ 



2669C(3) + 

12743tt 4 
1440 



6 

9095tt 6 
54432 



1826C(3) 2 4238C(5) 



+ W 



3 5 
96715C(3) 



53675 
24 



45 



+C\C F 



37684115 664325tt 2 
+ „„„ + 



5832 



1944 



18 

1265467vr 4 _ 
77760 



| 23tt 2 C(3) | 1616C(3) 2 



27 



18619tt 6 



272160 



1341553C(3) 355tt 2 C(3) 1136C(3) 2 2932C(5) 



52268375 



383660tt 2 
2187 



+ 



486 

152059tt 4 
19440 



27 
769tt 6 



102010C(3) 8149tt 4 
' + 3888 

17120104 



81 

+C A C F n f T 
416C(5) 



6561 

2 r7-i2 



+ 



+ C F njT 



5103 
343tt 2 C(3) 
27 

442961tt 2 

4374 ' 
2710864 
6561 ~ 



9 



+ C 2 F n f T 



9 

2732173 



1458 



13122 
45235tt 2 
486 



+ 



556C(5) 
45 

90148C(3) 
81 



1093tt 4 368tt 2 C(3) 



486 



27 



124tt 2 12784C(3) 



9 



243 



837T 4 

1215 



(17) 



i?(3) 
S 



r 3 



fin 

307057tt 2 
8748 



68590C(3) 77tt 2 C(3) 1766C(3) 2 20911C(5) 14474131 



+C>/T 



243 
8459vr 4 
+ 38880 
10021313 



+C F C A n f T 
~60 



6561 
155629 



108 
22523tt 6 
58320 
37868tt 2 
2187 



9 



45 



13122 



1508C(3) 437tt 4 _ 439?r 2 C(3) 6476C(5) 



27 1080 27 45 

64C(5)" 



+C F n f T 



+ 



243 9 
592C(3) 



41tt 2 23584C(3) 8vr 4 ,„ 2>/oN 
+ - — + 16vr 2 C(3) + 



81 



45 



9 



9 3 
2816C(3) 112vr 4 



320C(5) 



+ C F n)T 2 



42248 32tt 2 



9 



9 



135 



+ C A n)T 2 



2958218 152tt 2 47296C(3) 797tt 



6561 



+ 



81 



+ 



243 



+ 



1215 



,(18) 



where C F = (iV 2 - 1)/(2N C ), C A = N c , T = 1/2 and d abc d abc = (N 2 - 1)(N 2 - 4)/N c . 

We are confident that this technique can be applied to analytically evaluate master 
integrals appearing in various physical problems. 
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Appendix 

The coefficients of the dimensional recurrence relation for have the form 

8(d-4) 2 (d-3) 

° 7 ' 2 ~ ~5(d - 2)(d - l)(5d - 18)(5d - 16)(5d - 14)(5d - 12) ' 

(3d - 10) (483d 4 - 5996d 3 + 27684d 2 - 56272d + 42432) 
° 6 ' 3 ~ ~20(d - 2) 2 (d - l)(2d - 5)(5d - 18)(5d - 16)(5d - 14)(5d - 12) ' 
c 5j2 = _(d_3)[i5(d_4)(d-2) 2 (d- l)(2d-5)(3d- 10)(3d - 8)(5d - 18)]" 1 
x [(5d - 16)(5d - 14)(5d - 12)]" 1 [l2447d 7 - 256626d 6 + 2261972d 5 
-11052152d 4 + 32339200d 3 - 56684032d 2 + 55123200d - 22947840] , 
c 51 = - [60(d-4)(d-2) 2 (d-l)(2d-5)(3d-10)(5d-18)(5d-16)] _1 
x [(5d - 14) (5d - 12)] _1 [l8909d 7 - 384006d 6 + 3329804d 5 
-15982952d 4 + 45870976d 3 - 78731008d 2 + 74846208d - 30412800] , 

c 4 = _ [go(rf-3)(d-2) 2 (d- l)(3d- 10)(3d-8)(5d- 16)] _1 
x [(5d - 14) (5d - 12)] _1 [38619d 6 - 651987d 5 + 4575500d 4 
-17083884d 3 + 35791888d 2 - 39892032d + 18478080] . (19) 
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